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INTRODUCTION 


Waves  reaching  a  coastline  are  observed  to  break,  releasing  their 
energy  and  momentum  within  the  narrow  confines  of  the  surf zone.  This 
waveborne  momentum  flux  gives  rise  to  a  longshore  current,  which  con¬ 
tributes  to  the  longshore  sediment  transport  along  the  beach.  Under 
equilibrium  conditions,  the  longshore  sediment  transport  takes  place 
with  no  effect  on  the  existing  beach  bathymetry.  If  interrupted,  however, 
rapid  alteration  of  the  shoreline  can  occur,  which  may  jeopardize  near¬ 
shore  structures.  Recognizing  the  importance  of  these  processes  to  Navy 
operations,  the  Independent  Research  program  of  the  Naval  Civil  Engineering 
Laboratory  sponsored  a  study  to  develop  improved  models  to  predict  the 
magnitude  and  distribution  of  the  longshore  sediment  transport  rate  as  a 
function  of  incident  wave  and  beach  characteristics. 

Present  longshore  current  models  are  relatively  well-developed, 
being  based  on  a  momentum  balance  between  the  longshore  component  of 
waveborne  momentum  flux  and  the  bottom  and  lateral  shear  stresses  asso¬ 
ciated  with  the  longshore  current.  Early  longshore  current  models 
(Bowen,  1969;  Longuet-Higgins ,  1970)  assumed  near  normal  wave  angle 
incidence  and  weak  longshore  current  strength  in  order  to  simplify  their 
mathematical  solution.  Later  models,  such  as  those  by  Thornton  (1970), 
James  (1974),  Liu  and  Dalrymple  (1978),  and  Ostendorf  and  Madsen  (1979), 
assumed  arbitrary  sized  incident  wave  angles  and  longshore  current 
strengths . 

Present  longshore  sediment  transport  models  are  less  well-developed. 

In  general,  these  models  reflect  two  distinct  approaches:  an  energetics 
approach  and  a  traction  approach.  The  former  includes  models  by  Inman 
and  Bagnold  (1963),  Komar  (1971,  1977),  Thornton  (1973),  Bowen  (1981), 
and  Bailard  and  Inman  (1981),  while  the  latter  includes  models  by  Bijker 
(1971),  Swart  (1976),  and  Madsen  and  Grant  (1976).  In  both  cases,  these 
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approaches  are  based  on  adaptations  of  stream  flow  sediment  transport 
models.  The  energetics  approach,  which  is  used  in  this  report,  is  based 
on  the  stream  transport  model  developed  by  Bagnold  (1963,  1966). 

Bdgnold*s  energetics-based  sediment  transport  model  assumes  that 
the  sediment  is  transported  in  two  distinct  modes,  differing  by  the  way 
the  sediment  grains  are  supported.  Sediment  transported  as  bedload  is 
supported  by  the  bed  via  grain-to-grain  interactions,  while  sediment 
transported  as  suspended  load  is  supported  by  the  stream  fluid  via 
turbulent  diffusion.  In  both  modes,  energy  is  expended  by  the  stream  in 
transporting  the  sediment  load.  Bagnold,  comparing  the  stream  to  a 
machine,  defined  the  sediment  transport  efficiency  as  the  ratio  of  the 
rate  of  energy  expended  in  transporting  either  the  bedload  or  the  sus¬ 
pended  load  divided  by  the  total  rate  of  energy  production  of  the  stream. 
For  steady,  two-dimensional  stream  flow,  Bagnold  (1963)  developed 
the  following  total  load  sediment  transport  equation 

=  ^  ^ 

I  tan<t>  -  tanP  W  -  tanp 

u 

i  =  total  immersed  weight  sediment  transport  rate  (composed 
of  the  sum  of  the  bedload  transport  rate,  i^,  and  the 
suspended  load  transport  rate ,  ig) 

U)  =  rate  of  energy  production  of  the  stream 

u  =  mean  velocity  of  the  stream 

tanp  =  slope  of  the  stream  bed 

^  =  internal  angle  of  friction  of  the  sediment 

W  =  fall  velocity  of  the  sediment 

Cg  =  bedload  load  efficiency 

Eg  =  suspended  load  efficiency 

For  stream  flow  conditions,  Bagnold  (1966)  found  that  =  0.13  and 

eg  s  0.01. 


i 


where 
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For  oscillatory  flows,  such  as  those  found  in  the  surf zone,  Bagnold 
(1963)  developed  an  alternative  but  related  sediment  transport  equation. 
Bagnold  reasone.'  that  the  oscillatory  wave  motion  acts  to  move  the 
sediment  back  and  forth  in  an  amount  proportional  to  the  local  rate  of 
energy  dissipation.  Although  no  net  transport  results  from  this  motion, 
a  steady  current  of  arbitrary  strength,  when  superimposed  on  the  wave- 
induced  oscillatory  motion,  is  free  to  transport  the  sediment  in  the 
steady  current’s  direction.  This  conceptual  model  resulted  in  the 
following  sediment  transport  equation 


where 


=  K’  u) 


u 


(2) 


m 

u)  =  local  rate  of  energy  dissipation 

=  magnitude  of  the  oscillatory  water  velocity 

Uq  =  steady  current  in  the  0  direction 

=  time-averaged,  immersed  weight,  sediment  transport  rate  in 
the  6  direction 

K'  =  Bagnold* s  oscillatory  transport  rate  coefficient 


Equation  2  has  been  used  in  the  development  of  a  number  of  longshore 
transport  models,  including  those  by  Inman  and  Bagnold  (1963),  Komar 
(1971,  1977),  and  Thornton  (1973).  Inman  and  Bagnold’ s  model  is  currently 
the  most  widely  accepted  and  is  equivalent  in  form  to  the  equation 
recoomended  by  the  U.S,  Army  Corps  of  Engineers  (Anonymous,  1973). 

Inman  and  Bagnold* s  model  linearly  relates  the  total  spatially  integrated, 
immersed  weight,  longshore  transport  rate,  to  what  has  been  erroneously 

(Longuet-Higgins,  1972)  termed  the  longshore  component  of  wave  energy 
flux  per  unit  length  of  beach,  The  latter  is  defined  as 

=  (E  C  n)|^  cosa^  (3) 
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where  E  =  wave  energy  per  unit  length  of  crest 
C  =  wave  phase  velocity 
n  =  wave  energy  transport  function 
=  incident  breaker  angle 

b  =  subscript  that  denotes  the  value  at  the  breakpoint 

The  resulting  longshore  transport  equation,  termed  the  wave  power  equation, 
is  expressed  as 

=  K  Pjj,  (4) 

where  K  is  a  dimensionless  constant.  Komar  and  Inman  (1970)  found  that 
when  K  =  0.77,  Equation  4  best  fits  the  field  data  from  13  sand  tracer 
experiments  performed  at  El  Moreno  Beach,  Mexico,  and  Silver  Strand 
Beach,  California. 

Despite  the  general  acceptance  of  Equations  2  and  4,  their  develop¬ 
ment  is  only  loosely  based  on  Bagnold’s  stream  transport  model  (Equation  1). 
The  latter  suggests  that  an  attempt  should  be  made  to  differentiate 
between  the  bedload  and  suspended  load  transports  within  the  surfzone 
and,  in  addition,  the  local  bedslope  should  have  some  influence  on  the 
sediment  transport  rate.  Physical  reasoning  also  suggests  that  the 
time-averaging  procedure  implied  by  Equation  2  should  be  considered  more 
carefully.  Bailard  and  Inman  (1981)  and  Bowen  (1981)  have  attempted  to 
correct  some  of  these  difficulties.  Bailard  and  Inman  used  Bagnold*s 
stream-based  bedload  model  to  derive  a  general  form  bedload  equation 
that  is  valid  for  time-varying  flow  over  an  arbitrarily  oriented  sloping 
bottom.  Their  equation  reduces  to  Bagnold's  bedload  equation  for  steady, 
two-dimensional  stream  flow.  When  applied  to  the  surfzone,  Bailard  and 
Inman's  results  suggest  that  Equation  2  is  only  valid  for  weak  longshore 
currents  and  near  normal  wave  incidence  (when  suspended  sediment  transport 
is  neglected). 
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Bowen  (1981)  applied  Bagnold's  total  load  sediment  transport  equation 
to  the  problem  of  on-offshore  sediment  transport  on  beaches.  After 
generalizing  Bagnold's  total  load  Equation  1  for  the  special  case  of 
normally  incident  waves  and  steady  on-offshore  currents  (no  longshore 
currents),  Bowen  used  Stokes  second  order  wave  theory  and  Longuet-Higgins * 
(1953)  bottom  streaming  solutions  to  predict  the  equilibrium  beach 
profile  as  a  function  of  the  incident  wave  characteristics.  His  results 
support  observations  that  steep  beaches  are  coarse-grained  and  flat 
beaches  are  fine-grained.  Moreover,  small,  long-period  waves  cause  a 
beach  to  accrete  and  large,  short-period  waves  cause  a  beach  to  erode. 
Finally,  a  null-point  hypothesis  is  supported. 

At  the  present  time,  there  is  considerable  uncertainty  concerning 
the  relative  importance  of  bedload  versus  suspended  load  sediment 
transport  within  the  surf zone.  While  in  natural  streams  the  suspended 
sediment  transport  rate  generally  dominates  the  bedload  transport  rate, 
field  measurements  suggest  that  in  the  surfzone,  the  suspended  sediment 
transport  rate  accounts  for  less  than  20%  of  the  total  (Komar,  1978). 

These  measurements,  however,  implicitly  assume  that  all  sediment  trans¬ 
ported  within  approximately  10  cm  of  the  bottom  is  bedload.  In  stream 
flow,  however,  the  bedload  transport  is  generally  believed  to  take  place 
within  a  few  millimeters  of  the  stationary  bed.  As  a  consequence,  it 
would  seem  that  much  of  the  surfzone  sediment  transport  presently  con¬ 
sidered  to  be  bedload  may,  in  fact,  be  suspended  load.  This  hypothesis 
is  supported  by  Bowen’s  (in  press)  on-offshore  transport  model,  which 
suggests  that  the  relative  importance  of  the  two  sediment  transport 
rates  is  a  function  of  the  orbital  velocity  magnitude  divided  by  the 

sediment  fall  velocity,  u  /W. 

m 

The  purpose  of  this  report  is  to  present  a  total  load  (bedload  and 
suspended  load)  sediment  transport  model,  developed  for  time-varying 
flow  over  an  arbitrarily  sloping  planar  bed.  The  model  is  similar  to 
Bailard  and  Inman’s  (1981)  bedload  transport  model  in  that  it  predicts 
the  local,  nearbottom  sediment  transport  rate  as  a  function  of  the 
nearbottom  water  velocity  vector.  For  present  purposes,  nearbottom  is 
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considered  Lo  be  that  region  near  the  bottom  where  the  velocity  (iistri- 
bution  is  nearly  logarithmic  (see  the  Appendix).  The  model  is  similar 
to  Bowen^s  (1981)  model  in  that  Bagnold's  total  load  stream  transport 
model  (Equation  1)  is  used  as  a  working  hypothesis.  The  model  is 
specialized  for  idealized  surfzone  conditions  consisting  of  steady 
longshore  and  on-offshore  currents  in  addition  to  an  oscillatory  wave- 
induced  current  having  a  local  wave  angle.  The  model  is  found  to  predict 
the  local  equilibrium  beach  slope  as  well  as  the  local  longshore  sediment 
transport  rate.  In  addition,  the  model  is  combined  with  Ostendorf  and 
Madsen* s  (1979)  longshore  current  model  to  predict  the  magnitude  and 
distribution  of  the  longshore  sediment  transport  rate.  A  discussion  of 
the  model's  limitations  can  be  found  in  the  concluding  remarks. 


DEVELOPMENT  OF  A  TIME-VARYINO  TOTAL  LOAD  SEDIMENT  TRANSPORT  EQUATION 

Consider  the  time-varying  flow  of  water  with  vector  velocity  u^ 
moving  over  a  plane  sloping  bed  of  cohesionless  sediment.  Figure  1 
shows  the  x-axis  is  directed  upslope  and  the  y-axis  is  directed  cross¬ 
slope.  The  magnitude  of  the  slope  is  tanp,  and  the  velocity  u^  has  an 
arbitrary  orientation  relative  to  the  x-axis  but  is  coplanar  with  the 
x-y  plane.  The  shear  stress  on  the  bed,  T^,  is  assumed  to  be  described 
by 


where  p  =  water  density 

=  drag  coefficient  for  the  bed 
t  =  subscript  that  denotes  a  time-varying  quantity 


(5) 


Similarly,  the  local  rate  of  energy  dissipation,  is  assumed  to  be 
equal  to 


6 
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(6) 


U)^  =  p  Cf  |u^ 

Note  that  for  simplicity,  any  phase  difference  between  the  bed  stress 
and  the  nearbottom  fluid  velocity  is  neglected. 

Following  Bailard  and  Inman  (1981),  the  local  immersed  weight 
sediment  transport  rate,  i^,  is  assumed  to  be  a  time-varying  vector 
quantity  with  the  following  form 


where  and  Kg^  are  the  dimensionless  time-varying  vectors  associated 
with  the  bedload  and  suspended  load  transport  rates,  respectively. 
Assuming  the  bedload  transport  to  behave  as  a  granular-fluid  shear 
layer,  Bailard  and  Inman  (1981)  found  that 


where  =  bedload  efficiency 

D 

tan4>  =  internal  coefficient  of  friction  of  the  sand 

Equation  8  suggests  that  the  bedload  transport  rate  vector  is  composed 
of  two  components,  one  directed  parallel  to  the  instantaneous  velocity 
vector  and  the  other  directed  downslope. 

Bagnold*s  (1963,  1966)  conceptual  model  of  suspended  sediment 
transport  in  streams  is  relatively  simple.  He  reasoned  that  the  stream 
supports  the  suspended  sediment  via  turbulent  diffusion.  While  suspended, 
the  sediment  grains  move  downstream  with  nearly  the  local  fluid  velocity. 
At  the  same  time,  however,  they  are  also  falling  vertically  relative  to 
the  local  supporting  fluid  with  a  fall  velocity,  W.  In  equilibrium,  the 
center  of  mass  of  the  suspended  load  must  remain  at  a  constant  height 
above  the  stream  bed;  thus,  the  transport  process  extracts  energy  from 
the  stream  in  an  amount  equal  to  the  product  of  the  immersed  weight  of 
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the  suspended  sediment  and  the  sediment  fall  velocity.  Using  his  analogy 
of  the  stream  as  a  sediment  transporting  machine,  Bagnold  assumed  that 
the  power  expended  in  transporting  the  suspended  sediment  is  a  constant 
fraction,  £g,  of  the  total  power  produced  by  the  stream,  uj.  Extending 
this  conceptual  model  to  time-varying  flow  over  an  arbitrarily  sloping 
planar  bed,  it  can  be  shown  (see  the  Appendix)  that 


Similar  to  the  bedload  transport  rate,  the  suspended  sediment  transport 
rate  vector  is  seen  to  consist  of  two  components,  one  directed  parallel 
to  the  instantaneous  velocity  vector  and  the  other  directed  downslope. 
The  latter  is  associated  with  the  downslope  component  of  the  immersed 
weight  of  the  suspended  sediment. 

Combining  Equations  6,  7,  8,  and  9,  and  time-averaging  over  a  wave 
period  T,  where 


then  the  following  total  load  sediment  transport  equation  is  obtained 


P  c 


f  tan(t) 


tanP  ^ 
tan(t> 


/ 

3  ^ 

r 

“t 

u^>  -  —  tanP  < 

(11) 


Bagnold *s  stream-based  total  load  sediment  transport  equation  is  seen  to 
be  a  special  case  of  Equation  11,  since  both  equations  are  equal  when 

A 

tanp  <<  1  and  =  -u  i.  Equation  11  is  also  equivalent  to  Bowen's 
(1981)  initial  equation  for  on-offshore  sediment  transport  conditions 
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(Ut  = 


i),  except  for  the  extra  factor  £g  in  the  final  term.  This 


difference  is  of  particular  importance  in  relation  to  the  equilibrium 
beach  slope,  which  will  be  discussed  more  fully. 


LITTORAL  SEDIMENT  TRANSPORT 

Due  to  its  general  form,  Equation  11  is  not  particularly  useful 
unless  nearbottom  water  velocity  measurements  are  available.  These 
measurements  have  recently  become  possible  with  the  development  of 
robust  electromagnetic  current  meters  (Cunningham  et  al.,  1979).  For 
modeling  purposes,  however,  it  is  desirable  to  be  able  to  express 
Equation  11  in  terms  of  a  steady  current,  u,  oriented  at  an  angle  0  to 
the  beach  normal  and  an  oscillatory  wave-induced  velocity,  u,  oriented 
at  an  angle  a  to  the  beach  normal.  Figure  1  shows  the  instantaneous 
velocity  vector,  u^,  may  be  expressed  in  terms  of  its  x  and  y  components 

u^  =  (u  cosa  +  u  COS0)  i  +  (u  sinor  +  u  sin0)  j  (12) 

Observation  (Huntley,  1976)  has  shown  that  while  the  wave-induced  velocity 
is  by  definition  oscillatory  (<u>  =0),  it  is  not  symmetrically  distributed 
in  time.  Instead,  there  is  a  brief  high-velocity  onshore  flow  followed 
by  a  longer  but  less  intense  offshore  flow.  As  a  consequence,  the  third 
and  fifth  oscillatory  velocity  moments,  <u  >  and  <u  >,  are  nonzero. 
Following  Bailard  and  Inman  (1981)  or  Bowen  (1981),  this  asymmetry  may 
be  expressed  by  assuming  the  following  approximate  form  for  u 

u  =  u  cosat  •*’  u  ^  cos2at  +  .... 
m  m2 

where  u  >  u  etc.,  and  a  =  27t/T. 
m  mz 


(13) 


Equations  12  and  13  may  be  used  to  evaluate  the  time-averaged 
scalar  and  vector  quantities  in  Equation  11.  After  considerable  manipu¬ 
lation,  the  following  expressions  are  obtained  for  the  on-offshore, 

<i  >,  and  the  longshore,  <i  >,  sediment  transport  rates 
X  y 


<i  > 

X 


p  u 
^  f  m 


B 

tdn(|> 


3  (1  2  2 

4^,  cosa  +  6  +6  -r-  ^  cos  or  +  6 

^1  u  u  \  2  V 


+  6  sina  cosa  -  (u3)^ 

V  tan<^  ^  ^ 


u 

w-*^s 


4^2  cosa 


+  6  (u3)* 
u 


-  Cg  tanp  (u5)- 


(14) 


<i  > 

y 


p  c^  u 
^  f  m 


tan(t^ 


sina  +  6^  +  6  sin^a  +  6^ 

V  V  V  2  u 


+  6  sina  cosa 
u 


^  W“^S 


4^2  sina  +  6^(u3)* 


(15) 


where  the  relative  steady  current  strengths  6,  6^,  and  6^  are  defined  as 


u 

u 

m 

(16a) 

u 

T 

m 

COS0 

(16b) 

u 

IT 

m 

sin6 

(16c) 

The  velocity  moments  4^^  and  4>2  defined  as 

If/  =  <u^>/u^ 

m 


(17) 
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(18) 


^2  =  <|uP^>/u^ 

and  the  integrals  (u3)*  and  (u5)*  are  defined  as 


(19) 

(20) 

Bailard  and  Ininan  (1981)  obtained  approximate  analytic  expressions  for 
the  integral  (u3)*  for  the  special  cases  of  weak  and  strong  steady 
currents.  In  the  present  development,  however,  both  (u3)*  and  (u5)*  are 
left  unevaluated.  Later,  when  the  model  is  combined  with  Ostendorf  and 
Madsen's  (1979)  longshore  current  model  in  order  to  predict  the  distri¬ 
bution  of  the  longshore  sediment  transport  rate  as  a  function  of  the 
incident  wave  conditions,  these  integrals  are  evaluated  using  a  five-point 
Gaussian  quadrature  equation,  as  used  by  Ostendorf  and  Madsen,  i.e., 

(u3)*  =  0.118  [(u3)*(ot=2.99)  +  (u3)*(at=:0. 147) ] 

+  0.239  [(u3)*(at=2.42)  +  (u3)*(at=0.725)] 

+  0.284(u3)’^(at=1.57)  (21) 

Equation  21  can  be  shown  to  be  accurate  to  within  ±2%  for  6^2  and 
(6-a)  %  90  degrees.  Alternatively,  (u3)*  and  (u5)^  may  be  interpolated 
from  Figures  2  and  3. 

Equations  14  and  15  are  still  unsuitable  for  modeling  purposes 
because  of  the  presence  of  the  two  skewness  terms  and  Following 

Bailard  and  Inman  (1981),  however,  and  may  be  eliminated  from 
Equation  15  by  assuming  that  the  beach  is  in  local  on-offshore  equilibrium 
(<i^>  =  0).  Under  this  condition,  Equation  14  becomes 


/ 


(u3)*  =  V  /  +26  cos(0-a)  cosat  +  cos^at)^^^ 


dt 


(u5)*  ^  ^  J  +  26  cos(0-of)  cosat  +  cos^at)^^^ 


dt 
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4>i  cosa  = 


tan<^ 


(u3)*  -6^-6 
tani^  u  u 


1  +  cos^a  +  6^  j  -  6 

2  V  i  V 


sina  cosa 


IT  h  ^  t“3)* 


+  -^  tg  tanP  (u5)* 


(22) 


Combining  Equations  22  and  15 ,  one  obtains 


<i  > 

y 


p  C-  u 
^  f  m 


tan4» 


6-^  +  6^  6  +  + 
V  u  V  2 


iBfii  (u3)*  -  6^  6  -  6^  )  tana 

tan(|)  V  u  2  u  I 


u 


6^  (u3)*  -  6^(u3)*  tancf 


m  2 

+  “5" ^ »  tanP  (u5)*  tana 
VT  "" 


(23) 


Finally,  when  it  is  assumed  that  6^  <<  6^,  then  Equation  23  reduces  to 
the  desired  longshore  transport  equation 


<i  >  =  p  c-  u 

y  f  m 


tan<|> 


6^  +  ~  (u3)’^*  tana 

V  2  tan(t>  ^ 


u  ^0 

+  Cg  6^  (u3)*  +  Cg  tanP  (u5)*  tana 
w 


(24) 


Equation  24  can  be  compared  with  Bagnold's  oscillatory  transport 
Equation  2,  which  is  rewritten  for  longshore  transport  conditions  as 
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<i  >  =  p  K'  6 

y  3n  ^  f  m  V 


(25) 


Comparing  Equations  24  and  25 ,  both  equations  are  equal  if 


K’  =  ^ 


3n 


4  tan^ 


1+  6^  tana 

2  V  tan<^ 


4  W 


(u3)*  +  Cg  tanp  tana 

V 


(26) 


The  above  equation  suggests  that  instead  of  K'  being  constant,  it  is  a 
function  of  6^,  a,  tan^,  and  This  is  in  accordance  with  our 

earlier  hypothesis  based  on  Bagnold’s  stream  transport  (Equation  1). 
Mote  that  for  the  special  case  when  6^,  a,  tan^,  and  very 

small,  then  Bailard  and  Inman's  (1981)  result  is  recovered 


in 

8  tan^k 


(27) 


With  typical  surfzone  conditions  (tan^  =  0.10,  a  =  10  degrees,  =  0.21, 

and  Eg  =  0.025)  assumed,  Figure  4  shows  Bagnold's  transport  coefficient 

K'  as  a  function  of  the  dimensionless  longshore  current  strength,  5^, 

and  the  ratio,  u  /W.  The  coefficient  K*  is  seen  to  be  a  strong  function 
m 

of  both  6^  and  Only  for  weak  longshore  currents  and  small  values 

of  does  K'  approach  the  constant  value  expressed  by  Equation  27. 


EQUILIBRIUM  BEACH  SLOPES 

In  addition  to  predicting  the  longshore  sediment  transport  rate, 
the  present  model  can  be  used  to  predict  the  equilibrium  beach  slope  as 
well.  By  definition,  a  beach  slope  is  in  local  equilibrium  when  the 
divergence  of  the  time-averaged  sediment  transport  rate  vector  is  zero, 
i.e. , 
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5.<i  >  ^<i  > 

§— X-  +  i_L-  =  0 

dx  By  “ 


When  a  two-dimensional,  steady-state  beach  is  assumed,  then  ^  constant 

and  ^i^^  =  0*  H  the  beach  slope  tanp  is  small,  then  the  latter  condition 


applied  to  Equation  14  yields 


cosof  +  6^  +  6  6^  +  6  1’^  +  cos^a 

tan<J>  1  u  u  V  ^  \  2 


+  5  sina  cosa  +  cosa  +  6  (u3)^ 

V  W  S  ^2  u 


9 

-V  (u3)*  +-y4  (u5)* 
tan  <t)  Vr 


For  normal  wave  incidence,  negligible  longshore  currents,  and  weak 
on-offshore  currents,  Equation  29  becomes 


tanP  = 


tan<l)  (*’1  2  0  W  (*^2  *  3n  \ 


tan  <|) 


: _  .  16  % 

2.  15n  ^ 


Equation  30  is  closely  related  to  the  equilibrium  beach  slope 
equation  developed  by  Bowen  (1981).  The  two  equations  are  most  easily 
compared  if  the  skewness  factors,  and  and  the  relative  on-offshore 
current  strength,  6^,  are  estimated  from  Stoke* s  second-order  wave 
theory  and  Longuet-Higgins *  (1953)  bottom  drift  model,  respectively. 
Following  Bowen  (1981),  we  assume  that 


^  3  “m  -2 

u  =  u  cosat  +  T"7r“  sinh  kh  cos2at 
m  4  C 


and 


2 

u 

u  COS0  =  (32) 

where  C  =  wave  phase  velocity 
k  =  wave  number 
h  =  local  water  depth 


Combining  Equations  12,  17,  18,  31,  and  32,  the  skewnesses,  and 
are  found  to  be  equal  to 


.  7  U1 

*1^1  ~  1  r 


9  “m  .  . -2. . 
^  — sinh  kh 


and 


>1  12  %  .  . -2. ,  ^ 
‘*’2  =  C"  ^ 


(33) 


(34) 


and  the  dimensionless  on-offshore  current  strength,  6^,  is  equal  to 


5 

u 


(35) 


Combining  Equations  3C,  33,  34,  and  35,  the  equilibrium  beach  slope 
equation  then  becomes 


tanP 


tan6 


27n 

64 


^sinh  — 


tan^<|i 


5  w2"s 


sinh  kh  't’  4 


.  (36) 


The  principal  difference  between  Bowen's  equilibrium  beach  slope 

2 

equation  and  Equation  36  is  that  the  latter  has  an  in  the  denominator 
instead  of  an  e^.  As  a  result,  the  slope  effect  in  the  suspended  sediment 
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transport  rate  in  Equation  36  is  V^g  or  approximately  40  times  less 
than  in  Bowen’s  equation.  Predicted  equilibrium  beach  slopes  for  strong 
suspension  conditions  are  thus  significantly  greater  in  the  present 
model , 

Following  Bowen  (1981)  and  noting  that 


m  _  ak 
C  *  sin  kh 


(37a) 


u 

m  _  a  o 
W  ""  W  sinh  kh 


(37b) 


k 


tanh^^^  k  h 
o 


(37c) 


where  is  the  deepwater  wave  number,  then  the  modified  beach  slope, 
tanp/ak^,  may  be  plotted  as  a  function  of  the  parameters  k^h  and  aa/W 
(Figure  5),  The  resulting  figure  is  qualitatively  very  similar  to  the 
figure  obtained  by  Bowen.  Both  figures  show  that  the  beach  slope 
increases  with  decreasing  depth,  increasing  fall  velocity,  and  increasing 
wave  period.  All  three  features  are  generally  observed  in  nature.  A 
null  point  hypothesis  is  also  supported  in  that  sediment  coarser  than 
the  equilibrium  grain  size  moves  onshore,  while  sediment  finer  than  the 
equilibrium  moves  offshore.  The  biggest  discrepancy  between  the  present 
model  and  that  which  is  observed  in  nature  is  in  the  magnitude  of  the 
predicted  equilibrium  beach  slope  in  shallow  water.  The  present  model 
significantly  overpredicts  slopes  that  are  actually  observed  or  that  are 
predicted  by  Bowen’s  model.  The  evidence,  however,  is  inconclusive  in 
evaluating  either  model  because  Stoke ’s  second-order  wave  theory  is  not 
valid  in  shallow  water.  As  a  result,  it  is  best  to  consider  Equation  34 
or  Figure  5  as  indicative  only  of  trends  and  not  of  actual  magnitudes  of 
the  equilibrium  beach  slope. 

In  principle,  simultaneous  measurements  of  nearbottom  surf zone 
water  velocities  and  surfzone  beach  profiles  taken  over  a  considerable 
length  of  time  should  be  able  to  resolve  this  uncertainty.  A  preliminary 
analysis  of  field  data  from  the  Nearshore  Sediment  Transport  Study  at 
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Torrey  Pines  Beach  suggests  that  autosuspension  conditions  can  occur 
relatively  frequently  when  using  a  model  such  as  Bowen’s,  which  is 
compatible  with  Bagnold's  equations.  This,  of  course,  leads  to  momentary 
infinite  offshore  suspended  sediment  transport  rates,  which  seems  unten¬ 
able.  The  present  model  avoids  this  problem  for  reasons  discussed  in 
the  Appendix.  A  more  comprehensive  study  of  this  question  is  needed, 
however,  before  any  conclusions  can  be  made. 

The  above  discussion  pertained  only  to  normally  incident  waves  with 
negligible  longshore  currents.  Equation  29,  however,  permits  an  evalua¬ 
tion  of  the  relative  effects  of  wave  angle  incidence  and  longshore 
current  strength  on  the  local  equilibrium  beach  slope.  Due  to  the 
complexity  of  Equation  29,  it  is  necessary  to  consider  nominal  surfzone 
conditions  where  =  0.21,  ”  0.025,  =  0.10,  =  0.10,  and 

tan())  =  0.63.  Figure  6  shows  tanp  plotted  as  a  function  of  the  relative 
longshore  current  strength,  6^,  and  the  local  wave  angle,  a.  It  is 
evident  that  for  these  conditions,  tanp  is  a  strong  function  of  the 
longshore  current  strength,  decreasing  with  increasing  6^  for  6^  >  0.2. 
The  equilibrium  beach  slope  is  a  somewhat  weaker  function  of  the  local 
wave  angle,  a  increasing  slightly  with  increasing  a  for  6^  >  0.04, 

The  sensitivity  of  the  equilibrium  beach  slope  to  the  longshore 
current  strength  suggests  a  mechanism  for  the  formation  of  a  breakpoint 
bar  and  trough.  Field  and  laboratory  observations  of  longshore  currents 
suggest  that  the  longshore  current  strength  increases  with  distance  from 
the  shoreline  until  reaching  a  maximum  near  the  breakpoint.  It  subse¬ 
quently  decreases  with  additional  distance  offshore.  Figure  6  suggests 
that  the  effect  of  the  longshore  current  on  the  equilibrium  beach  profile 
would  be  to  flatten  the  profile  in  the  area  of  the  maximum  longshore 
current  (i.e.,  just  inside  the  breakpoint).  This  flattening  would,  in 
effect,  create  a  bar  and  trough  type  of  feature,  with  the  bar  located 
seaward  of  the  current  maximum  and  the  trough  located  shoreward  of  the 
maximum. 


SURFZONE  SEDIMENT  TRANSPORT  DISTRIBUTION  AND  INTEGRATED  TOTAL 


Equation  15  predicts  the  local  total  longshore  transport  rate  as  a 
function  of  the  local  wave-induced  velocity  magnitude,  the  local 

wave  angle,  Of,  and  the  local  longshore  current,  6^.  In  order  to  predict 
the  distribution  of  the  longshore  sediment  transport  rate  across  the 
surfzone  as  a  function  of  the  incident  wave  conditions,  it  is  necessary 
to  incorporate  a  longshore  current  model.  A  review  of  existing  longshore 
current  models  suggests  that  Ostendorf  and  Madsen's  (1979)  model  is  most 
suitable  for  the  present  study.  The  important  features  of  their  model 
are  that  finite  wave  height,  longshore  current  strength,  and  incident 
wave  angle  effects  are  incorporated,  as  are  the  effects  of  lateral 
momentum  exchange. 

Ostendorf  and  Madsen's  (1979)  longshore  current  model  uses  a  simi¬ 
larity  form  of  solution,  where  the  fundamental  solution  is  approximately 
equal  to  Longuet-Higgins '  (1970)  linear  solution,  and  correction  factors 
are  included  to  account  for  the  effects  of  finite  longshore  current 
strength,  incident  wave  angle,  and  breaker  height.  The  model  has  the 
following  form 

V  =  X  n^  cosa^ 

where  ^  linear  longshore  current  distribution 

X  =  finite  current  strength  factor 

n^  =  factor  equivalent  to  the  wave  energy  transport  function 
at  the  breakpoint 

cosa^  -  finite  incident  wave  angle  factor 

The  dimensionless  longshore  current  distribution,  v*,  has  the  following 
form 

/  ^  ■  s  \  s 

V*  =  Cj  I  y-Tc  )  +  Cj  X*  (?  *  0.4,  (39) 
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V* 

V* 

V* 

where 


C 


C 


C 


1 


2 


3 


0,26  X*  -  0.71  X*  In  x* 
0.26 


(1  -  2.5  P) 


(P  t  0.4,  x*>l) 


(P  =  0.4,  x*Sl) 


(P  =  0.4,  x*>l) 


(40) 


The  dimensionless  distance  from  the  shoreline,  x*  (see  Figure  1), 
is  defined  as 


x*>  = 


(41) 


where  Xg  and  Xg  are  the  positions  of  the  shoreline  and  breakpoint, 
respectively.  The  lateral  mixing  parameter,  P,  is  defined  as 


P 


71  r  tanA 
2  c. 


(42) 


where  F  is  a  lateral  mixing  coefficient  and  tanA  is  equal  to  the  average 
beach  slope,  tanP,  modified  to  account  for  the  wave  setup,  i.e., 


tanA 


tanP  (1  +  |- 


(43) 


Based  on  laboratory  velocity  profiles,  Ostendorf  and  Madsen  found  that 
r  ^  0.13. 
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In  deriving  Equation  38,  it  is  assumed  that  the  local  wave-*height* 
to-water-depth  ratio,  y  =  H/h,  has  the  following  distribution 


Y 


(x*>l) 


(44) 


In  addition,  the  characteristic  longshore  current  magnitude,  v^,  is 
defined  as 


V 

c 


6  u  - 
c  mb 


(45) 


where  the  magnitude  of  the  oscillatory  wave  velocity  at  the  breakpoint, 
u^j^,  is  obtained  from  shallow  water  wave  theory,  i.e.. 


“mb  = 


(46) 


and  the  relative  longshore  current  strength,  6^,  is  defined  as 


5 

c 


5n  tanA  sinof. 


8  c. 


(47) 


Ostendorf  and  Madsen  present  an  iterative  procedure  for  estimating  the 
breakpoint  values  of  H,  h,  and  n  as  a  function  of  the  incident  wave 
conditions.  In  addition,  the  finite  current  parameter.  A,  is  obtained 
iteratively  from  the  parameters  n^,  a^,  6^,  and  P, 

Combining  the  longshore  current  Equation  38  and  the  longshore 
transport  Equation  24,  and  noting  that 
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(48) 


m 


mb  Y 


X 


and 


% 


(49) 


then  the  longshore  transport  rate  distribution  becomes 


<i  >  =  p  c-  u  , 

y  f  mb 


tan^ 


COSOb  U]  V*  X* 


+  6^  cos^Ob  V*^  +  (u3)* 


tana 


+  -5^  € 

W 


6^  \  Ob  cosoib  (u3)*  (—  I  v*x* 


,3/2 


+  ^  tanp 


(u5)* 


y  V  „.5/2 


X*  '  tana 


(50) 


Equation  50  can  be  normalized  by  dividing  by  the  bedload  transport  rate 
at  the  breakpoint  under  weak  current,  near  normal  wave  incidence  condi¬ 
tions  with  no  lateral  mixing,  i^,  where 


1  B  3 

i  =  T  p  c  -  r — r  5  ,  cosa, 

c  2  f  tan<^  D  c  mb  b 


(51) 


Following  this  procedure,  Equation  50  becomes 
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* 

<i  > 
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<i  > 
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=  A  I  X )  V*  X*  +  2  6^  A^  cos^a  v*^ 

CD  D 


tan^ 


X  T — A  (u3)*(-^l  tana 

6c  coscr^^  tan^  I  y  ' 


jLr  .-p/2 


“mb 

^  -vT 


2  A  (u3)*  -i_  v*x* 


2  .2 

“mb 

+  — ^  -- —  tan<t)  tanp 

vr 


6  n,  cosa.  \  y, 
c  b  b  \  'b 


)  (u5)*  tana 


(52) 


Because  of  its  complexity,  Equation  52  is  difficult  to  generalize 

in  a  simple  manner.  Nevertheless,  assuming  nominal  surfzone  conditions 

(tanp  =  0.1,  \  =  1.0,  n^  =  1.0,  =  10  degrees,  6^  =  0.5,  P  =  0.2, 

£g  =  0.21,  and  8g  =  0.025),  the  distribution  of  can  be  plotted  as  a 

function  of  the  ratio  (Figure  7).  In  Figure  7,  several  features 

can  be  seen.  First,  the  profiles  are  a  strong  function  of  the  parameter 

u^^/W.  Increasing  values  of  strongly  increases  the  total  load 

sediment  transport  rate.  Another  feature  of  the  profiles  is  the  slope 

discontinuity  at  the  breakpoint  (x*  =  1).  This  feature  is  a  curiosity 

associated  with  the  slope  discontinuity  in  the  distribution  of  y  at  this 

point.  In  reality,  this  discontinuity  would  not  exist.  A  final  feature 

of  the  profiles  is  that  the  relative  location  of  the  sediment  transport 

rate  maximum  is  seaward  of  the  maximum  longshore  current. 

The  ratio  of  the  normalized  suspended  sediment  transport  rate 

divided  by  the  normalized  bedload  transport  rate,  i^^/i  is  also  a 

S  B 

function  of  the  distance  from  the  beach.  Figure  8  shows  a  plot  of  the 

distribution  of  the  ratio  i  */i  *  as  a  function  of  the  ratio  u  , /W  for 

o  B  mb 

the  same  surfzone  conditions  as  those  assumed  for  Figure  7.  As  might  be 
expected,  the  ratio  increases  with  distance  offshore  until  the  breakpoint, 
and  then  begins  to  decrease.  Another  anticipated  result  is  that  the 
ratio  increases  throughout  the  nearshore  with  increasing  values  of 

“mb/'^- 
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A  more  easily  measured  quantity  than  the  longshore  transport  rate 
distribution  is  the  spatially  integrated  total  longshore  transport  rate, 
defined  as 

X 

s 

I.  =  /  <i  >  d 

£  y  y  X 


tanA 


<i*>dx* 

y 


(53) 


Substituting  Equations  50,  51,  and  52  into  Equation  53  and  noting  that 


"b  8 


P  g 


(54) 


and 


(C  n). 


% 


y,  mb 


(55) 


then  Equation  53  becomes 


Y  £  00 

“TT  f  (E  C  n),  sina,  cosO, 

16  tan()>  J  y  b  b  I 

0 


(56) 


Comparing  Equation  56  with  the  wave  power  Equation  4,  the  wave  power 
coefficient  K  is  seen  to  have  the  following  form 


K 


5n  \  h 
16  tan(t> 


dx* 


(57) 


Instead  of  being  constant,  Equation  57  suggests  that  K  is  a  complex 

function  of  the  incident  wave  conditions,  the  beach  slope,  and  the 

sediment  size.  For  the  special  case  of  negligible  lateral  mixing,  weak 

longshore  current  strength,  near  normal  breaker  angle,  small  breaker 

height,  and  small  water-velocity-to-fall-velocity  ratio,  K  becomes  a 

constant,  K  ,  where 
*  o’ 
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5n  "b 

LO  ^  ^ 


(58) 


Again  the  normalized  wave  power  coefficient,  K*,  is  a  complex 
function  of  surfzone  parameters.  To  illustrate  its  behavior,  we  assume 
nominal  surfzone  conditions  where  tanA  =  0.1,  n^  =  1.0,  A  =  1.0,  =  1.0, 

P  =  0.2,  Cg  =  0.21,  and  Cg  =  0.025.  Figure  9  shows  K“^  plotted  as  a 
function  of  the  relative  longshore  current  strength,  6^,  the  breaker 
angle,  Of^^,  and  the  ratio  small  values  of  is  seen 

to  be  nearly  independent  of  6^  or  As  the  ratio  increases, 

however,  K*  also  increases,  becoming  a  function  of  6^  and  a^.  Finally, 

K*  decreases  with  increasing  longshore  current  strength  and  increases 
with  increasing  breaker  angle.  These  results  will  vary,  however,  for 
different  input  conditions. 

The  relative  overall  fraction  of  suspended  sediment  transport 
versus  bedload  transport,  also  be  plotted  as  a  function  of 

6^,  Of^,  and  Assuming  the  same  surfzone  conditions,  Figure  10 

shows  that  the  relative  transport  ratio  behaves  exactly  as  the  normalized 
wave  power  coefficient,  a  result  to  be  expected  because  for  large  values 
of  most  of  the  sediment  transport  is  by  suspension. 


MODEL  CALIBRATION 

Laboratory  and  field  data  available  for  calibrating  and  evaluating 
the  longshore  transport  model  are  relatively  limited.  Greer  and  Madsen 
(1978)  reviewed  existing  field  data  and  concluded  that  Komar  and  Inman's 
(1970)  sand  tracer  studies  at  El  Moreno  and  Silver  Strand  Beaches  were 
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the  most  reliable.  Very  recently,  a  series  of  large-scale  field  experi¬ 
ments,  termed  the  Nearshore  Sediment  Transport  Study,  have  been  conducted 
at  two  Southern  California  beaches.  Details  of  the  studies  may  be  found 
in  Cable  (1979,  1981);  however,  the  data  from  these  studies  are  only  now 
becoming  available  for  analysis  and  thus  will  not  be  considered  in  the 
present  study , 

Laboratory  data  are  also  very  limited.  In  the  present  case,  measure 
ments  must  include  the  incident  wave  conditions,  the  mean  longshore 
current,  the  equilibrium  beach  slope,  and  the  total  longshore  sediment 
transport  rate.  Following  the  study  by  Ostendorf  and  Madsen  (1979), 
selected  data  from  Saville  (1950)  and  Shay  and  Johnson  (1951)  will  be 
used  in  the  present  study.  A  summary  of  the  field  and  laboratory  data 
used  in  the  present  study  may  be  found  in  Table  1. 

The  numerical  calibration  procedure  is  as  follows.  In  reference  to 
Equations  50,  52,  and  56,  the  present  longshore  transport  model  has 
three  unknown  parameters:  the  drag  coefficient,  c^;  the  bedload  effi¬ 
ciency,  and  the  suspended  load  efficiency,  Using  Ostendorf  and 

Madsen’s  longshore  current  model,  the  drag  coefficient  is  estimated  from 
the  measured  mean  longshore  current.  In  the  process,  estimates  for  E^, 
^c^  ^Iso  obtained.  These  parameters 

completely  define  the  longshore  current  distribution.  Next,  the  distri¬ 
bution  of  <i*>  is  calculated  as  a  function  of  €«  and  using  Equation  52 
y  Do 

Integrating  this  distribution  via  Equation  57,  the  estimated  wave  power 

coefficient,  K  .  ,  is  calculated  as  a  function  of  e-,  and  yielding 
es  t  D  o 
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^  f  2  A  (u3)*  ty  V*  x*3/2 
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(u5)^'  tana 


dx*  (63) 


Estimates  of  K^,  K^,  and  for  all  data  sets  are  included  in  Table  1. 

With  use  of  these  estimates,  a  nonlinear  least-square  procedure  (Draper 
and  Smith,  1966)  was  used  to  estimate  the  values  of  the  bedload  efficiency, 
Cg,  and  the  suspended  load  efficiency,  £g,  from  the  measured  values  of 
the  wave  power  coefficient,  •  This  procedure  requires  constructing 

a  two-dimensional  contour  plot  of  the  mean  square  error,  S(eg,  £g) * 
defined  as 


=  =  rh  \  '>'obs-''est)' 

j=l 

where  n  is  the  total  number  of  data  pairs. 

While  following  the  above  procedures,  it  was  found  that  the  data 
from  El  Moreno  Beach  appeared  to  be  anomolous.  When  compared  with  the 
other  data,  the  measured  values  of  K  appeared  to  be  roughly  three  times 
too  large.  As  a  consequence,  it  was  decided  to  calibrate  the  model 
using  only  the  field  data  from  Silver  Strand  Beach  and  the  labor.itory 
data.  The  resulting  least-square  contour  plot  is  shown  in  I igure  11 
The  least-square  estimates,  £g  and  e^,  were  found  to  be  equal  to  0.?' 
and  0.025,  respectively.  Figure  11  also  shows  the  95%  confidence  contour, 
suggesting  that  the  95%  bounds  on  and  £-  are  0  <  C-.  <  0.44  and 

Do  n 

0.016  <  £g  <  0.031.  Figure  12  shows  a  plot  of  the  relative  errors, 
^obs^^est*  ^  function  of  the  estimated  wave  power  coefficient. 

Clearly  the  relative  error  of  the  calibration  data  is  approximately 
equally  distributed  about  one,  suggesting  that  £g  and  are  unbiased 
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Figure  12  also  shows  the  relative  errors  of  the  El  Moreno  Beach 
data.  These  errors  show  that  the  measured  wave  power  coefficients  are 
seriously  underestimated.  One  explanation  for  this  behavior  is  that  the 
transport  processes  at  El  Moreno  Beach  may  have  been  significantly 
different  than  those  at  Silver  Strand  Beach  and  in  the  laboratory.  In 
particular,  the  steep  beach  face  and  small  wave  heights  produced  waves 
which  broke  on  the  beach  face  itself.  It  would  be  expected  that  swash 
zone  transport  processes  would  predominate  under  these  conditions. 
Observations  have  shown  that  the  swash  zone  is  a  region  of  abnormally 
high  suspended  sediment  concentrations  and  transport,  and  as  a  result, 
it  appears  plausible  that  the  present  model  would  tend  to  underestimate 
transport  in  this  region. 

The  least  square  estimates,  and  are  encouraging  because  they 
are  similar  in  size  to  the  values  estimated  by  Bagnold  for  stream  flow. 
Bagnold  (1966)  analyzed  data  from  laboratory  and  field  measurements  of 
sediment  transport  rates  in  channel  flow  and  estimated  that  =  0.13 

D 

and  Cg  =  0.01.  The  fact  that  the  present  estimates  are  roughly  twice 
these  values  is  not  improbable.  In  the  case  of  the  bedload  efficiency, 
most  of  the  transport  in  the  surfzone  is  suspension,  and  the  total 
transport  rate  is  relatively  insensitive  to  the  estimated  value  of  the 
bedload  efficiency.  This  is  reflected  in  the  relatively  large  error 
bounds  on  Bagnold* s  stream-based  estimate  clearly  falls  within  the 

allowable  error.  In  the  case  of  the  suspended  sediment  efficiency, 

Bagnold* s  estimate  does  not  fall  within  the  bounds  of  the  present  estimate. 
This  may  reflect  the  different  turbulence  structure  of  the  surfzone 
relative  to  the  stream  flow.  In  fact,  the  surfzone  turbulence  simulation 
model  developed  by  Johns  (1980)  suggests  that  the  turbulence  energy 
density  near  the  bed  is  dominated  by  the  turbulence  energy  production  at 
the  face  of  the  bore  as  opposed  to  that  generated  by  the  shear  at  the 
bed.  The  latter  is  clearly  the  dominant  process  in  stream  flow.  As  a 
result,  it  is  not  surprising  that  the  suspended  sediment  efficiency  is 
enhanced  in  the  surfzone. 
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CONCLUDING  REMARKS 


The  present  sediment  transport  model  incorporates  a  wide  range  of 
simplifying  assumptions.  While  these  assumptions  are  necessary  due  to 
the  complexity  of  actual  surfzone  processes,  it  is  believed  that  they 
have  some  physical  basis. 

A  number  of  assumptions  used  in  the  present  model  share  a  common 
basis  with  many  longshore  current  models.  These  assumptions  include: 
the  use  of  linear  shallow  water  wave  theory,  linear  breaker  height 
variation  across  the  surfzone,  and  a  bed  stress  formulation  which  is 
Quadratic  in  velocity  and  has  a  zero  phase  lag.  Some  of  these  assumptions 
are  better  justified  than  others.  Guza  and  Thornton  (1980)  have  shown 
that  linear  wave  theory  does  a  surprisingly  good  job  of  predicting 
oscillatory  velocity  magnitudes  from  local  water  depth  measurements. 
Laboratory  measurements  (Bowen  et  al.,  1968)  have  shown  that  for  normal 
incident  waves,  a  linear  breaker  height  variation  across  the  surfzone  is 
reasonably  accurate.  The  quadratic  stress  formulation,  however,  appears 
more  open  to  question.  Laboratory  and  theoretical  studies  of  purely 
oscillatory  turbulent  boundary  layers  have  shown  that  the  nearbottom 
water  velocity  lags  the  bed  shear  stress  in  phase  by  approximately 
30  degrees  (Fisher  et  al.,  1974;  Jonsson  and  Carleson,  1976;  Jonsson, 

1980;  Smith,  1977;  Grant  and  Madsen,  1979).  The  applicability  of  these 
studies  inside  the  surfzone  is  unknown.  The  bore  propagation  model 
developed  by  Johns  (1980),  however,  appears  to  be  more  realistic  of 
surfzone  conditions.  This  model  suggests  that  a  quadratic  shear  stress 
model  with  a  constant  drag  coefficient  underestimates  the  magnitude  of 
peak  bed  stresses  which  occur  under  the  face  of  the  bore.  Despite  these 
problems,  the  quadratic  shear  stress  model  is  utilized  in  most  longshore 
current  models  because  it  is  well-suited  to  a  vertically  integrated 
formulation. 

The  present  energetics-based  sediment  transport  model  also  uses 
vertically  integrated  equations.  As  a  consequence,  no  details  of  the 
vertical  distribution  of  suspended  sediment  flux  are  resolved,  A  further 
consequence  of  using  vertically  integrated  equations  is  that  the  sediment 
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transports  are  assumed  to  respond  to  the  nearbottom  water  velocity  in  an 
instantaneous,  quasi-steady  manner.  This  assumption  is  probably  valid 
for  bedload  transport  (except  for  a  phase  lag)  because  the  thickness  of 
the  bedload  layer  is  small  and  can  respond  quickly  to  the  instantaneous 
bed  shear  stress.  The  suspended  sediment  transport,  however,  is  dis¬ 
tributed  over  a  greater  layer  thickness,  typically  on  the  order  of 
several  centimeters.  Assuming  a  layer  thickness,  J0,  and  a  sediment  fall 
velocity,  W,  the  characteristic  time  constant  for  the  sediment  transport 
layer  is  £/W,  which  is  typically  on  the  order  of  1  to  2  seconds.  For 
most  planar  beaches  with  incident  wave  periods  of  roughly  8  seconds,  it 
appears  that  the  quasi-steady  assumption  is  reasonable. 

For  rippled  bottoms,  the  situation  is  more  complicated.  Inman  and 
Bowen  (1963)  found  that  in  the  laboratory,  simple  phase-independent 
sediment  transport  models  are  inadequate  for  combined  oscillatory  and 
steady  flow  over  a  rippled  bottom.  They  found  such  flows  to  exhibit  a 
complex  phase-dependent  behavior  associated  with  vortex-induced  suspension 
occurring  over  the  ripple  crests.  As  a  consequence,  the  present  model 
is  probably  most  useful  in  describing  sediment  transport  inside  the 
surfzone,  where  plane  bed  conditions  usually  prevail.  Limiting  the 
model  to  this  region  also  minimizes  the  effects  of  neglecting  the 
threshold  stress  condition  for  sediment  transport  due  to  the  stronger 
steady  currents  found  there. 

Other  uncertainties  of  the  present  model  include  the  use  of  constant 
bedload  and  suspended  load  efficiency  factors.  Although  constant  values 
have  been  found  to  be  adequate  in  steady  stream  flow  (Bagnold,  1966), 
their  adequacy  in  time- varying  flow  is  unknown. 

Despite  the  above  limitations,  the  present  model  does  exhibit 
certain  behavior  observed  in  nature.  This  is  most  evident  in  its  pre¬ 
dictions  of  the  equilibrium  beach  slope.  The  present  study  substantiates 
Bowen’s  (1981)  results;  coarse-grained  beaches  are  steeper  than  fine¬ 
grained  beaches  and  short-period,  large-amplitude  waves  are  erosive.  It 
is  also  encouraging  that  the  present  model  produces  results  similar  to 
the  wave  power  equation.  To  a  large  degree,  however,  this  is  the  result 
of  both  approaches  being  based  on  energy  considerations.  Finally,  the 


fact  that  the  bedload  and  suspended  load  efficiency  factors  estimated 
from  surfzone  data  are  similar  in  size  to  the  values  estimated  from 
stream  flow  lends  further  credence  to  the  model. 

In  closing,  several  comments  can  be  made  concerning  improvements  to 
the  model.  It  is  clear  from  the  above  discussion  that  many  of  the 
details  of  the  sediment  transport  processes  are  hidden  by  the  use  of 
vertically  integrated  equations.  In  part,  this  is  justified  by  the 
resulting  mathematical  simplicity  and  by  the  lack  of  adequate  data  to 
validate  more  sophisticated  boundary  layer  models.  The  latter,  however, 
is  changing  due  to  the  advent  of  more  sophisticated  sensors.  Clearly, 
if  significant  progress  is  to  be  made,  models  which  predict  the  details 
of  the  boundary  layer  flow  field  and  sediment  concentration  distribution 
will  have  to  be  developed  for  the  surfzone.  If  simple  predictive  equations 
are  desired,  it  is  always  best  to  develop  simplifying  assumptions  based 
on  a  knowledge  of  the  complete  solution. 


REFERENCES 

Anonymous  (1973).  Shore  protection  manual,  U.S.  Army  Corps  of  Engineers, 
Coastal  Engineering  Research  Center.  1973. 

Bagnold,  R.  A.  (1963).  "Mechanics  of  marine  sedimentation,"  in  The  Sea; 
Ideas  and  Observations,  volume  3.  New  York,  N.Y. ,  Interscience  Publishers, 
1963,  pp  507-528. 

Bagnold,  R.  A.  (1966).  An  approach  to  the  sediment  transport  problem 
from  general  physics,  U.S.  Geological  Survey,  Professional  Paper  No. 

422-1.  Washington,  D.C.,  1966. 

Bailard,  J.  A.,  and  D.  L.  Inman  (1981).  "An  energetics  bedload  transport 
model  for  a  plane  sloping  beach;  local  transport,"  Journal  of  Geophysical 
Research,  vol  86,  no.  3,  1981,  pp  2035-2043. 


30 


Bijker,  E.  W.  (1971).  **Longshore  transport  computations/^  Proc.  ASCE 
Journal  of  Waterways,  Harbors  and  Coastal  Engineering,  WW4,  Nov  1971. 

Bowen,  A.  J.  (1969).  ’’The  generation  of  longshore  currents  on  a  plane 
beach,”  Journal  of  Marine  Research,  volume  27,  no.  2,  1969,  pp  206-215. 

Bowen,  A.  J.  (1981).  ”Simple  models  of  nearshore  sedimentation;  beach 
profiles  and  longshore  bais,”  in  Proceedings  of  the  Conference  on  Coastline 
of  Canada,  Halifax,  1978,  Geological  Survey  of  Canada. 

Bowen,  A.  J. ,  D.  L.  Inman,  and  V.  P.  Simmons  (1968).  ”Wave  set-down  and 
set-up,”  Journal  of  Geophysical  Research,  vol  73,  no.  8,  1968,  pp  2569 -2577. 

Cunningham,  P.  M. ,  R.  T.  Guza,  and  R.  L.  Lowe  (1979).  "Dynamics  cali¬ 
bration  of  electromagnetic  flow  meters,”  I.E.E.E.  Oceans,  vol  79,  1979, 
pp  298-301. 

Draper,  N.  R. ,  and  H.  Smith  (1966).  Applied  regression  analysis.  New 
York,  N.Y.,  John  Wiley,  1966. 

Fisher,  J.,  B.  Johnson,  and  M.  E.  McCormick  (1974).  "Time  dependent 
shear  stress  beneath  a  shoaling  wave,”  in  Proceedings  of  the  International 
Symposium  on  Ocean  Wave  Measurements  and  Analysis,  New  Orleans,  1974, 
vol  1.  New  York,  American  Society  of  Civil  Engineers,  1974,  pp  799-816. 

Gable,  C.  G.  (1979).  Report  on  data  from  the  Nearshore  Sediment  Transport 
Study  experiment  at  Torrey  Pines  Beach,  California,  November  -  December 

1978,  Institute  of  Marine  Research,  IMR  Ref.  No.  79-8.  La  Jolla,  Calif., 

1979. 


Gable,  C.  G.  (1981).  Report  on  data  from  the  Nearshore  Sediment  Transport 
Study  experiment  at  Ledbetter  Beach,  Santa  Barbara,  CA,  January  -  February 
1980,  Institute  of  Marine  Research,  IMR  Ref.  No.  80-5.  La  Jolla,  Calif., 


Grant,  W.  D.  ,  and  0.  S.  Madsen  (1979).  ’’Combined  wave  and  current 
interaction  with  a  rough  bottom,”  Journal  of  Geophysical  Research, 
vol  84,  no.  C4,  1979,  pp  1797-1808. 

Greer,  M.  N. ,  and  0.  S.  Madsen  (1978).  ’’Longshore  sediment  transport 
data;  A  review,”  in  Proceedings  of  the  16th  Conference  on  Coastal 
Engineering.  American  Society  of  Civil  Engineers,  1978. 

Guza,  R.  T.,  and  E.  B.  Thornton  (1980).  "Local  and  shoaled  comparisons 
of  sea  surface  elevations,  pressures  and  velocities,”  Journal  of 
Geophysical  Research,  vol  85,  no.  C3,  1980,  pp  1524-1530. 

Huntley,  D.  A.  (1976).  "Lateral  and  bottom  forces  on  longshore  currents,” 
in  Proceedings  of  the  15th  Conference  on  Coastal  Engineering,  Honolulu, 
Hawaii,  1976.  American  Society  of  Civil  Engineers,  pp  645-659. 

Inman,  D.  L. ,  and  R.  A.  Bagnold  (1963).  "Littoral  processes,”  in  The 
Sea;  Ideas  and  Observations,  volume  3.  New  York,  N.Y.,  Interscience 
Publishers,  1963,  pp  529-533. 

Inman,  D.  L. ,  and  A.  J.  Bowen  (1963).  "Flume  experiments  on  sand  transport 
by  waves  and  currents,”  in  Proceedings  of  the  Eighth  Conference  on 
Coastal  Engineering.  American  Society  of  Civil  Engineers,  1963,  pp  137-150. 

James,  I.  D.  (1974),  "A  nonlinear  theory  of  longshore  currents,”  Estuarine 
and  Coastal  Marine  Science,  vol  2,  1974,  pp  235-299. 

Johns,  B.  (1980).  "The  modeling  of  the  approach  of  bores  to  a  shoreline,” 
Coastal  Engineering,  vol  3,  no.  3,  1980,  pp  207-220. 

Jonsson,  I.  G.  (1980).  ”A  new  approach  to  oscillatory  rough  turbulent 
boundary  layers,”  Ocean  Engineering,  vol  7,  no.  1,  1980,  pp  109-152. 


32 


JoQSson,  I.  G.,  and  N.  A.  Carleson  (1976).  ''Experimental  and  theoretical 
investigations  in  an  oscillatory  rough  turbulent  boundary  layer,"  Journal 
of  Hydraulic  Research,  vol  lA,  1976,  pp  45-60, 

Komar,  P.  D.  (1971),  "The  mechanics  of  sand  transport  on  beaches," 

Journal  of  Geophysical  Research,  vol  76,  no.  3,  1971,  pp  713-721. 

Komar,  P.  D.  (1977).  "Beach  sand  transport:  Distribution  and  total 
drift,"  Proceedings  of  ASCE  Journal  of  Waterways,  Harbors  and  Coastal 
Engineering,  WWA,  May  1977,  pp  225-239. 

Komar,  P.  D.  (1978).  "Relative  quantities  of  suspension  versus  bedload 
transport  on  beaches,"  Journal  of  Sedimentary  Petroleum,  vol  48,  no.  3, 
1978,  pp  921-932. 

Komar,  P.  D. ,  and  D,  L.  Inman  (1970).  "Longshore  sand  transport  on 
beaches,"  Journal  of  Geophysical  Research,  vol  75,  no.  30,  1970,  pp  5914- 
5927. 

Lui,  P.  L-F. ,  and  R.  Dalrymple  (1978).  "Bottom  frictional  stresses  and 
longshore  currents  due  to  waves  with  large  angles  of  incidence,"  Journal 
of  Marine  Research,  vol  37,  no,  2,  1978,  pp  357-375. 

Longuet-Higgins ,  M.  S.  (1953).  "Mass  transport  in  water  waves," 
Philosophical  Transactions  of  the  Royal  Society,  London,  Series  A, 
vol  245,  1953,  pp  535-581. 

Longuet-Higgins,  M.  S.  (1970).  "Longshore  currents  generated  by  obliquely 
incident  sea  waves,  2,"  Journal  of  Geophysical  Research,  vol  75,  no.  33, 
1970,  pp  6790-6801. 


Longuet-Higgins,  M.  S.  (1972).  "Recent  progress  in  the  study  of  longshore 
currents,"  in  Waves  on  Beaches  and  Resulting  Sediment  Transport,  R.  E. 
Meyer,  Editor.  New  York,  N.Y. ,  Academic  Press,  1972,  p.  462. 


Madsen,  0.  S.,  and  W.  D.  Grant  (1976).  Sediment  transport  in  the  coastal 
environment,  Massachusetts  Institute  of  Technology,  Department  of  Civil 
Engineering,  Ralph  M.  Parsons  Laboratory,  Report  No.  241.  1976,  p.  105. 

Ostendorf,  D.  W. ,  and  0.  S.  Madsen  (1979).  An  analysis  of  longshore 
currents  and  associated  sediment  transport  in  the  surfzone,  Massachusetts 
Institute  of  Technology,  Department  of  Civil  Engineering,  Ralph  M. 

Parsons  Laboratory.  Cambridge,  Mass,,  1979,  p.  169. 

Saville,  T.  (1950).  "Model  study  of  sand  transport  along  an  infinitely 
long,  straight  beach,"  Transactions,  American  Geophysical  Union,  vol  31, 
no.  4,  1950,  pp  555-565. 

Shay,  E.  A.,  and  J.  W.  Johnson  (1951).  Model  studies  on  the  movement  of 
sand  transported  by  wave  action  along  a  straight  beach,  University  of 
California,  Department  of  Engineering,  Institute  of  Engineering  Research, 
Series  14,  Issue  7.  1951,  p.  58. 

Smith,  J.  D.  (1977).  "Modeling  of  sediment  transport  on  continental 
shelves,"  in  The  Sea:  Ideas  and  Observations,  volume  6,  New  York, 

N.y. ,  Interscience  Publishers,  1977,  pp  539-577. 

Smith,  J-  D.,  and  S.  R.  McLean  (1977).  "Spatially  averaged  flow  over  a 
wavy  bottom,"  Journal  of  Geophysical  Research,  vol  82,  no.  12,  pp  1735-1746. 

Swart,  D.  H.  (1976).  "Predictive  equations  regarding  coastal  transports," 
in  Proceedings  of  the  15th  Conference  on  Coastal  Engineering.  American 
Society  of  Civil  Engineers,  1976. 

Thornton,  E.  B.  (1970).  "Variation  of  longshore  current  across  the 
surfzone,"  in  Proceedings  of  the  12th  Conference  on  Coastal  Engineering. 
American  Society  of  Civil  Engineers,  1970,  pp  291-308. 

Thornton,  E.  B.  (1973).  "Distribution  of  sediment  transport  across  the 
surfzone,"  in  Proceedings  of  the  13th  Conference  on  t.onstal  Engineering. 
American  Society  of  Civil  Engineers,  1973,  pp  1049-'1068. 


34 


Table  1.  Data  Relating  Longshore  Sediment  Transport  to  Wave  and  Beach  Characteristics 
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Figure  6.  Relationship  between  the  equilibrium  beach  slope,  tan/5,  and  the 
dimensionless  longshore  current  strength,  6^,  for  different  local 
wave  angles,  Q. 
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Figure  7.  Dimensionless  surf^o^e  total  load  sediment  transport  and  velocity 
profiles  for  different  values  of  the  velocity  ratio, 


Figure  8.  Surfzonc  distribution  of  the  sedimciU  transport  ratio  i^/i^,  for 
different  values  of  the  velocity  ratio 
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Figure  10.  Reiationship  between  the  ratio  of  the  sediment  transport  ratio,  I5/I 
and  the  relative  longshore  current  strength,  6^,  for  different  values 
of  the  velocity  ratio 
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igure  12.  Comparison  between  the  estimated  and  observed  wave  power  coefficients 
for  several  laboratory  and  field  studies. 


Schematic  diagram  used  in  developing  a  general  form  suspended 
sediment  transport  model. 


Appendix 


A  GENERAL  DERIVATION  OF  BAGNOLD’S 
SUSPENDED  SEDIMENT  TRANSPORT  EQUATION 


For  simplicity,  we  first  consider  steady  two-dimensional  stream 
flow  over  a  planar  bed  with  slope  tanp.  The  water  flows  with  velocity 
u(z)  over  a  depth  h.  The  x-coordinate  is  directed  downslope  and  the 
z-coordinate  is  directed  perpendicular  and  upward  from  the  bed.  The 
origin  is  at  the  base  of  the  bed. 

From  geometry,  the  rate  of  energy  production  of  the  stream, 
is  equal  to 


h 

u)  ^  -fog  sinB  u*i  dz  (Al) 

stream  J  ®  ^ 

0 

where  p  is  the  apparent  density  of  the  sediment-ladened  water.  This 
a 

may  be  expressed  as 

P.  =  (P„  '  p)N  +  p  (A2) 

where  =  density  of  the  sediment  grains 

N  =  local  volume  concentration  of  solids 
p  =  density  of  sediment-free  water 

Assuming  that  the  suspended  sediment  moves  downstream  with  approximately 
the  local  water  velocity,  u,  then  Equations  Al  and  A2  may  be  combined  to 
yield 
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0  ^  =  (p  *  p)  g  sin3  /n  u*i  dz  + 

stream  '^s  # 


where  U)  is  the  sediment-free  stream  power. 

Following  Bagnold,  we  define  an  immersed  v^eight  suspended  sediment 
transport  rate  per  unit  area  bed,  ig,  as 


lo  =  (p^  -  p)  8  cosp  /n  u  dz 

i>  s  Q 


Bagnold  (1963,  1966)  assumed  that  the  rate  of  energy  dissipation  per 
unit  area  of  bed  associated  with  this  transport  of  sediment, 
equal  to  the  sediment  load  times  its  fall  velocity,  i.e., 


‘"sed  =  (Ps  ■  8 


where  W  is  the  sediment  fall  velocity.  Combining  Equations  A4  avj  AS, 
we  find 


where 


/ »: 
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6agnold*s  critical  hypothesis  is  that  the  rate  of  energy  dissipation 
associated  with  the  suspended  sediment  transport,  related  to 

the  total  rate  of  energy  production  of  the  stream,  through  a 

constant  suspended  load  efficiency  factor,  £g,  so  that 


U) 


sed 


stream 


(A8) 


Combining  Equations  A3,  A4,  A6,  and  A8,  we  obtain  an  equation  for  the 
magnitude  of  the  suspended  sediment  transport  rate 


(A9) 


It  should  be  noted  that  Equation  A9  does  not  reduce  to  Bagnold*s  stream- 
based  transport  model  except  for  negligible  slopes.  To  see  this,  we 
rearrange  Equation  A9  to  become 


eg  U. 


-  Cg  tan^^ 


(AlO) 


By  inspection,  Equation  AlO  differs  from  Bagnold's  stream  transport 
Equation  1  by  the  inclusion  of  the  factor  £g  in  the  denominator.  Con¬ 
ceptually  this  difference  arises  because  Bagnold  assumed  that  the  stream 
power  contribution  from  the  suspended  sediment  load  contributes  directly 
to  the  suspended  sediment  transport  rate  instead  of  through  an  efficiency 
factor  as  in  Equation  A9.  If  Bagnold* s  hypothesis  is  correct,  then 
Equation  A9  should  read 


'f 


=  ^  (tan^  +  tg  u)) 


(All) 
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This  subtle  difference  between  Equations  A9  and  All  is  negligible 
for  the  sma 1 1  he<ls lopes  norma  1 ly  in  stream  flow  condi tions .  It  becomes 
extremely  important,  however,  with  respect  to  equilibrium  beach  slopes 
as  seen  in  the  main  body  of  the  paper.  This  is  because  Equation  All  (or 
Equation  1)  predicts  autosuspension  conditions  (infinite  downslope 
transport  rates)  when  [uj  exceeds  W/tanp.  Interestingly,  this  condition 
can  be  exceeded  momentarily  on  natural  sand  beaches.  For  autosuspension 
to  occur  in  the  present  model,  however,  |u[  must  exceed  W/(eg  tanp) . 

Since  ^  ^0,  this  condition  is  much  less  apt  to  be  met. 

Although  Equation  A9  was  developed  for  two-dimensional  stream  flow, 
it  seems  plausible  tiiat  it  can  be  modified  to  include  the  more  general 
situation  of  a  steady  flow  over  a  planar  sloping  bed  having  an  arbitrary 
orientation  relative  to  the  direction  of  fluid  flow.  Under  these  more 
general  conditions,  ihe  direction  of  the  sediment  transport  rate  need 
not  be  coiinear  with  the  direction  of  fluid  flow.  In  addition  to  moving 
with  the  nearbottom  water  velocity,  the  downslope  component  of  suspended 
sediment  load  can  be  expected  to  produce  a  downslope-directed  sediment 
transport . 

In  the  general  case  of  suspended  sediment  transport,  sedimenc  can 
be  expected  to  be  distributed  throughout  the  turbulent  boundary  layer 
region.  In  the  case  of  fully  developed  stream  flow,  this  is  roughly  the 
total  water  depth,  but  for  oscillatory  flow  the  boundary  layer  thickness 
is  much  less.  Smith  (1977)  and  Grant  and  Madsen  (1979)  show  that  the 
turbulent  boundary  layer  thickness  under  waves  is  approximately  (u*)^/a, 
where  is  the  magnitude  of  the  oscillatory  shear  velocity  and  a  is 

the  angular  frequency  of  the  waves.  For  typical  nearshore  wave  conditions, 
the  boundary  layer  thickness  is  on  the  order  of  several  centimeters. 

Within  the  boundary  layer  region,  the  sediment  concentration  increases 
exponentially  towards  the  bed.  Consequently,  most  of  the  sediment  is 
found  relatively  near  the  bed.  For  reasons  that  will  become  evident,  it 
is  matheinatically  convenient  to  assume  that  most  of  the  sediment  is 
transported  near  the  bed,  within  the  logarithmic  or  *'law  of  the  wall*' 
portion  of  the  boundary  layer,  where  the  shear  stress  is  constant.  For 
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steady  two-dimensional  stream  flow,  the  logarithmic  boundary  layer  can 
be  assumed  to  be  the  lower  one-tenth  of  the  stream  depth  (Smith  and 
McLean,  1977).  For  oscillatory  turbulent  boundary  layers,  Smith  (1977) 
concludes  that  the  logarithmic  layer  thickness  is  0.02  times  the  wave 
boundary  layer,  while  Grant  and  Madsen  (1979)  assume  that  it  is  the  full 
boundary  layer  thickness  as  defined  by  The  latter  will  be 

assumed  in  the  present  development. 

Consider  a  steady  current,  u,  flowing  over  a  planar  bed  with  slope 
tanp  (Figure  13).  The  current  flows  parallel  to  the  bed  and  at  an  angle 
0  relative  to  the  downslope-directed  x-axis.  The  y-axis  is  directed 
across  slope  and  the  z-axis  is  directed  perpendicular  and  upward  from 
the  bed.  The  origin  is  assumed  to  be  at  the  top  of  the  bedload  layer, 
which  is  nominally  a  few  grain  diameters  above  the  stationary  bed.  As 
discussed,  the  suspended  sediment  transport  is  assumed  to  occur  relatively 
near  the  bed,  as  compared  with  the  total  water  depth.  Although  in 
practice  the  local  transport  rate  will  be  colinear  with  the  local  shear 
stress,  it  is  assumed  for  simplicity  that  the  direction  of  the  suspended 
sediment  transport  rate  is  co-linear  with  the  shear  stress  at  the  base 
of  the  suspended  sediment  layer. 

Considering  now  the  limited  nearbottom  region  of  suspended  sediment 
transport  and  assuming  that  steady,  nonaccelerating  conditions  exist, 
the  momentum  equation  is 


(Pg  -  P)N  I  +  =  0 


(A12) 


where  ?  is  the  sediment-ladened  fluid  stress  tensor,  which  has  been 
corrected  for  the  sediment-free  fluid  pressure,  pgz.  Boundary  conditions 
for  the  flow  are  as  follows.  At  the  top  of  the  suspension  layer,  the 
normal  stress,  is  zero.  At  the  base  of  the  suspension  layer,  the 

current  exerts  a  shear  stress,  T^,  in  the  6  direction.  Assuming,  however, 
that  the  suspension  layer  is  sufficiently  small,  the  current-related 
shear  stress  can  be  assumed  to  be  applied  to  the  top  of  the  suspension 
layer.  This  is  equivalent  to  assuming  that  the  current-related  shear 
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stress  is  constant  throughout  the  suspension  layer,  as  is  consistent 
with  a  logarithmic  velocity  profile.  Assuming  that  the  logarithmic 
boundary  layer  thickness  is  £,  the  boundary  conditions  may  be  expressed 
as 
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where  is  the  downslope  component  of  stress  due  to  the  immersed  weight 
of  the  suspended  sediment  and  is  equal  to 


Tg  =  (p^  -  p)  g  sinP  j  H  dz 


(A19) 


Assuming  that  <<  T^,  which  is  consistent  with  the  assumption  of  a 
small  bedslope^  then 


f  -  =  T  +  T  COS0 

z=0  o  g 


(A20) 


?  T  -  T 

—  =  (cos0  +  sin^6)i  +  (sin6  -  sin6  cos6)  j  (A21) 

1^|  2=0  'o  ^0 


Noting 9  however,  that 


— .  =  cosB  i  +  sinS  j 


T  ^  r>  ^  ^ 

,-r]  ~  r^l  (sin^e  i  -  sinO  cosB  j) 

|T|  2=0  |u|  o 


(A22) 


As  previously  meationed,  the  suspended  sediment  transport  rate  is 
assumed  to  be  directed  parallel  to  ^^-0’ 


(A23) 
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Combiniag  Equations  A9,  A22»  and  A23,  we  obtain  to  o(e|) 


t 


W 


^  s 

u)  u  +  — r  tanp  u) 


cos6 


u 


(sin^8  i  -  sin6  cos0  j) 


where  it  has  been  assumed  that  u  =  |  u 
Assuming  next  that 


(A24) 


T  = 
o 


and 


U)  = 


(A25) 


(A26) 


where  is  the 


drag  coefficient  for  the  bed, 


then  Equation  A24  becomes 


•t 


P  -vT 


+  p  c^  tanP 


(A27) 


Assuming  now  that  the  nearbottom  velocity  vector  is  a  time-varying 
quantity,  u  ,  and  that  the  suspended  sediment  transport  rate  responds  to 
the  instantaneous  velocity  in  a  quasi-steady  manner,  a  vector,  may 

be  defined  as 


(A28) 
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Thus, 


+  Eg  tanP 


(A29) 


Finally,  if  the  bed  slopes  in  the  negative  i-direction,  Equation  A29 
becomes  equal  to  Equation  9  in  the  main  body  of  text. 
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NOMENCLATURE 


s  -  s 


X 


-t 


<i  > 

X 


<i  > 

y 


Wave  amplitude  (£) 

Phase  velocity  (il/t) 

Modified  Longuet-Higgins  constants 
Bed  drag  coefficient 
Wave  energy  (m/t^) 

Gravity  (£/t^) 

Wave  height  (£) 

Water  depth  (JH) 

Local  sediment  transport  rate  (m/t^) 

Local  sediment  transport  rate  vector  (ro/t^) 

Unit  vector  in  onshore  direction 

Spatially  integrated  longshore  bedload  sediment  transport 
rate 

Local  bedload  sediment  transport  rate  (m/t^) 

Local  bedload  sediment  transport  rate  vector  (m/t^) 

Characteristic  local  longshore  sediment  transport  rate 
(m/t^) 

Spatially  integrated,  total  load,  longshore  sediment 
transport  rate  (m-£/t^) 

Spatially  integrated  longshore  suspended  load  sediment 
transport  rate  (m-£/t^) 

Local  suspended  sediment  transport  rate  (m/t^) 

Local  suspended  sediment  transport  rate  vector  (m/t^) 

Time-averaged  local  on-offshore  sediment  transport  rate 
(m/t^) 

Time-averaged  local  longshore  sediment  transport  rate 
(m/t^) 
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<i*> 

y 


y 

j 

K 

K* 


N 


n 


% 

P 


P 


£ 


S(eB»  ^5) 


Normalized  time*- averaged  local  longshore  sediment 
transport  rate 

Local  sediment  transport  rate  in  6  direction  (m/t^) 

Unit  vector  in  longshore  direction 
Wave  power  coefficient 

BagnoId*s  oscillatory  transport  rate  coefficient 
Normalized  wave  power  coefficient 
Wave  power  coefficient  integrals 
Vector  bedload  wave  power  coefficient 
Estimated  wave  power  coefficient 

Wave  power  coefficient  for  negligible  lateral  mixing, 
longshore  current  strength,  breaker  angle,  wave  height, 
and  orbital-velocity-to-fall“velocity  ratio 

Observed  wave  power  coefficient 

Vector  suspended  load  wave  power  coefficient 

Wave  number  (l/£) 

Deep  water  wave  number  (l/£) 

Sediment  layer  thickness  (£) 

Volume  concentration  of  solids 

Wave  energy  transport  function 

Finite  breaker  height  factor 

Lateral  mixing  parameter 

Longshore  component  of  wave  energy  flux  (m*£/t^) 

Mean  square  error  of  wave  power  coefficient  estimates 
Wave  period  (t) 

Sediment-ladened  fluid  stress  tensor  (m/t^) 

Downslope  component  of  immersed  sediment  weight  (m/£*t^) 
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zz 


u 

(u3)* 

(u5)* 

“e 

V 

V* 

V 

C 

W 

X 

X* 

*s 

y 

Z 

GKEEK 

O 

“b 

tanp 

r 


Shear  stress  at  the  base  of  the  suspension  layer  (m/£*t^) 
Normal  stress  at  top  of  the  suspension  layer  (m/£-t^) 

Time  (t) 

Mean  velocity  magnitude  (£/t) 

Oscillatory  surfzone  velocity  (£/t) 

Total  vector  velocity  (£/t) 

Velocity  magnitude  integral 
Velocity  magnitude  integral 
Oscillatory  velocity  magnitude  (£/t) 

Oscillatory  shear  velocity  magnitude  (£/t) 

Mean  water  velocity  in  0  direction  (£/t) 

Longshore  current  magnitude  (£/t) 

Dimensionless  longshore  current  distribution 
Characteristic  longshore  current  magnitude  (£/t) 

Sediment  full  velocity  (£/t) 

Shore  normal  coordinate  (£) 

Dimensionless  surfzone  position 
Shoreline  position  (£) 

Breaker  position  (£) 

Longshore  coordinate  (£) 

Vertical  coordinate  (£) 

ALPHABET 

Wave  angle  (deg) 

Breaker  angle  (deg) 

Bed  slope 

Lateral  mixing  coefficient 
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fM 


A 

6 

6 


0 

\ 

P 

P, 


‘l'2 

u> 

U) 

sed 

stream 


Ratio  of  wave  height  to  water  depth 
Bed  angle  modified  for  wave  setup  (deg) 

Dimensionless  steady  current  strength 
Relative  longshore  current  strength 
Dimensionless  onshore  current  strength 
Dimensionless  longshore  current  strength 
Bedload  efficiency 
Suspended  load  efficiency 
Steady  current  angle  (deg) 

Finite  longshore  current  factor 
Fluid  density  (m/£^) 

Apparent  density  of  suspension  (m/£^) 

Sediment  grain  density 
Wave  frequency  (1/t) 

Bed  shear  stress  vector  (m/£-t^) 

Internal  angle  of  friction  (deg) 

Velocity  skewness  factor 
Velocity  skewness  factor 

Local  rate  of  energy  dissipation  production  (m/t^) 

Rate  of  energy  dissipation  associated  with  the  suspended 
sediment  transport  (m/t^) 

Rate  of  energy  production  of  a  stream  (m/t^) 


OTHER  SYMBOLS 


<  >  Time  average  quantity 

I  I  Absolute  value 

Vector  quantity 
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SUBSCRIPTS 


Break  point  quantity 
Time  varying  quantity 


b 
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COMOCEANSYSPAC  SCE.  Pearl  Harbor  HI 
C’OMSUBDEVGRUONE  Operations  Offr.  San  Diego.  CA 
NAVSl.^RFLANT  Norfolk.  VA 
NAVSIJRFPAC  Code  N-4,  Coronado 

COMOPTEVFOR  CMDR.  Norfolk.  VA:  Code  705.  San  Diego,  CA 

DEFFUELSITPCEN  DFSC'-OWE  (Grafton)  Alexandria.  VA;  DFSC^OWE  (Term  Engrng)  Alexandria.  VA 

DLSIE  Army  Logistics  Mgt  Center.  Fort  Lee,  VA 

DNA  STTL.  Washington  DC 

OTIC  Defense  Technical  Info  Clr'Alexandria.  VA 

DTNSRDC  Anna  Lab  (Code  119)  Annapolis  MD;  Anna  Lab  (Code  L5('i8)  Annapolis  MD 
FMFLANT  CEC  Offr.  Norfolk  VA 
FMFPAC  CG(FEO)  Camp  Smith.  HI 

OS  A  Assist  Comm  Dcs  &  C'nst  (FAIA)  D  R  Dibner  Washington,  DC 
LIBRARY  OF  C’CJNGRESS  Washington.  DC  {Sciences  &  Tech  Div) 

MARINE  CORPS  BASE  )st  For  Serv  Supp  Gru  (C'SS'5)  Camp  Pendleton  CA;  PWO.  C  amp  S.  D.  Butler, 
Kawa.saki  Japan 

MCAS  Facil.  Engr.  Div.  Cherry  Pom!  NC;  Facs  Main!  Dept  -  Operations  Div.  Cherry  Point 
MCDEC  M&L  Div  Quantico  VA;  NSAP  REP,  Ouantico  VA 
MCRD  PWO.  San  Diego  Ca 

MILITARY  SEALIFT  COMMAND  Washington  DC 
NAF  PWO.  Atsugi  Japan 

NARF  Eituipmeni  Engineering  Division  (Code  610(K)).  Pensacola,  FL 

NAS  Dir.  Util.  Div..  Bermuda;  Lead.  Chief.  Petty  Offr.  PW/Self  Help  Div.  Bccville  TX;  PW  (,L  Maguire), 
Corpus  Chrisii  TX;  PWD  -  Engr  Div.  Oak  Harbor.  WA;  PWD  Maint.  Div.,  New  Orleans.  Belle  Chasse  LA: 
PWn.  Code  1821H  (PfanKuch)  Miramar.  SD  CA;  PWO  Belle  Chasse,  LA;  PWO  Key  West  FL;  PWO. 
Cilcnvicw  IL;  ROICC  Key  West  FL.  SCE  Norfolk,  VA:  Shore  Facil.  Ofr  Norfolk.  VA 
NATL  BliREAU  OF  STANDARDS  Kovacs.  Washington.  D  C. 

NATL  RESEARCH  COUNCIL  Naval  Studies  Board.  Washington  DC 
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NA\  AC’^I  PWO.  lAMidon  VK 
NAVArROSPRECiMElX'l  N  St'E.  Pensacola  FL 
N  A\  AIRDEA'CFN  Code  81.C  Warminsler  PA 
NAXC'HAPCfRC  CO  Wdliamslnirf’  VA 

NANCOASrSYSCFN  Code  172  (C  H  Koesv )  Panama  Cits  M 

N AX'C’OASI SVSTC.' FR  C’O.  Panama  C  it\  F  I  .  Code  (.1  (^uirk)  IFmama  CH\  )i  ;  (  ode  (J  Ntitllcmani 
Panama  Citv.  Kl  ,  labrarv  Panama  City.  FI. 

NAVC'OMMARFAMSI  RSl  A  SCE.  Wahiawa  FFl:  SCI  I  mi  1  Naples  IiaK 

N.AVC’OMMS FA  Code  41)1  Nea  Makri.  Greece.  PWD-F  NCt  l)i\  I^ir  Washineion.  IH  .  PWO,  l-.\niouih. 
.Australia 

NAVC'ONS  FRACF-N  Curriculum  Inslr  Sids  Ollr.  (iullporl  MS 
NAVEOFRAPRODEVC'EN  Technical  I.ihrars.  FVnsacola.  FI 
NAVEl.EXSYSC'OM  Code  F'ME  i:4-{d.  Washington.  IX  .  PMF  04  oi:.  Wjsh  DC 
NAVEODFAC'  CYide  Indian  Head  MD 
\A\  FAC'  PWO.  C  enterville  Bch.  Ferndale  C  A 

NA\  FAC'ENCiC'O.M  Cinle  ItTF  (Essoglou)  .Alexandria.  \'A.  C  ode  n4'  .Mex.tiuln.i.  '  \  {  ynW  044  Alcx.indii.i, 
YA;  Code  0451  (P  W’  Brewer)  .Alexandria.  \’.i;  C'udc  04^^  iD  i  Alcx.iiulii.i,  \  \.  (  lule  o4>X  , 

.Alexandria.  YA:  CTide  IMA  I  Alexandria.  \’ A;  C  tide  i)o.  Alexandn,}  \  A  t  .ule  ooO  .Mcx.indria.  \  .A:  (  ode 
111.'  (M.  C  arr)  .Alexandria.  \’A:  C'ode  111.'  (1.  Stevens)  Alexandiia.  \  \.  \li'i!tNOii  X  ,»p  (  ,iu*)ine  K 
NAVF.ACENGCOM  -  CHES  DIV.  Code  40.S  Wash.  IX  ;  Code  40^  (1)  SchccsGci  Wadiingu<n  IX  .  C  ode 
FFO-lC'  Washington  l^C':  FPO-I  Washington.  DC',  FFXJ-IFi.A.^  Washingi»m  (H 
NA\TACE.NGCT)M  -  1  ANI  DIV.  Eur  BR  Deputy  Dir.  Naples  ltal>.  RDIM  FC)  lo:  Noitolk  \  . A 
NAVFAC'ENCiCOM  -  NORTFI  DIV.  (Borctskv)  Philadelphia.  PA;  CO.  C\h1c  oop  Phiiadelplii.i  PA.  Code  lols. 

RDTAELC).  Philadelphia  PA:  ROICC.  Contracts,  C  rane  IN 
NAVFACENGC'O.M  -  PAC'  DIV.  CODE  IWP  PEART  HARBOR  HI;  C  ode  :oil  Pearl  llarhor.  HI:  C  ode  40:. 

RDT&E.  Pearl  Harbor  FFI:  Commander,  Pearl  Harbor.  HI 
NAVFACENGCOM  -  SOUTH  DIV.  Code  90.  RDT&ELC4,  Charleston  SC 

NAVFACENGCOM  -  WEST  DIV.  Code  04B  San  Bruno.  C’A;  O9P20  San  Bruno.  C  A.  RD!A:EI.O  Code  2011 
San  Bruno.  CA 

NAVFACENGCOM  CONTRACF  Eng  Div  dir.  Southwest  Pae.  Manila.  PI:  OICC.  Southwest  Pac.  Manila.  PI: 

OICC^ROFCC.  Balboa  Panama  Canal;  ROICC,  KeOavik,  Iceland:  ROICC.  NAS,  Corpus  Christi,  TX 
NAVFORCARFB  Commander  (N42).  Puerto  Rico 
NAVMAG  SCE.  Subic  Bay,  R.P. 

NAVOCEANO  Library  Bay  St.  Louis,  MS 

NAVOCEANSYSCEN  Code  41.  San  Diego,  CA;  Code  4475  Bayside  Library,  San  Diego,  CA;  Code  447,'IB 
(Tech  Ub)  San  Diego.  CA;  Code  52  (H.  Talkinglon)  San  Dicgt>  C'A:  Code  5204  (.1  Slaehiw).  San  Diego. 
CA;  Code  5214  (H.  Wheeler).  San  Diego  CA;  Code  5221  (R. Jones)  San  Diego  C'a;  Code  5311  (Bachman) 
San  Diego,  CA;  Hawaii  Lab  (R  Yumori)  Kailua.  HI;  Hi  Lab  Tech  Lib  Kailua  HI 
NAVPGSCOL  C.  Morers  Monterey  CA;  Code  bIWL  (O.  Wilson)  Monterey  C'A;  E  Fluirnltm.  Monierex  CA 
NAVPHIBASE  CO.  ACB  2  Norfolk,  VA;  COMNAVBEACHGRCi  FWO  Norfolk  \  A:  C  ode  S-3F.  Noiiolk  VA 
Dir.  Amphil'*.  Warfare  Brd  Staff.  Norfolk,  VA;  Harbor  C'learance  Unit  Fwiv  Little  C  reek.  \',\:  SC'li. 
Coronado.  San  Diego  C'A 
NAVREGMEDCEN  SCE  (D.  Kaye);  SCE.  (iuam 
NAVSCOLCECOFF  C35  Port  Mueneme,  CA 
NAVSeSOL  PWO.  Athens  GA 

NAVSEASYSCOM  C'ode  SEA  OOC'  Washington.  DC':  SEA  (ME  (L  Kess)  Washington.  DC’ 

NAVSEC’GRUAC'7  PW'O.  Adak  AK 

NAVSHIPYD  Bremerton.  WA  (Carr  Inlet  Acoustic  Range);  Code  202.5  (F.ibrary)  Puget  Sound.  iXe  inert  on 
WA;  C'ode  .380.  Portsmouth.  VA;  C'ode  440  Portsmouth  NFL  C  ode  440.  I’ugel  Sound.  Bremerton  W'A:  leeh 
Library,  Vallejo.  C'A 

NAVSTA  Dir  Engr  Div,  PWD.  Mayport  FL;  PWD  (l.TJCi.P.M.  Motolenich).  FYieito  Rico;  PWO.  Maxport  FI.. 

Utilities  Engr  Off  Rota  Spain 
NAVTECHTRACEN  SCE.  Pensacola  FL 

NAVWPNC'EN  Code  26.^6  (W.  Bonner).  China  Lake  C'A:  Code  3803  China  l.akc,  C'A 
NAVWPN.STA  PW  Office  Yorktown.  VA 

NAVWPNSTA  PWD  -  Maint.  Control  Div.,  Concord.  CA;  PWO.  C  harleston.  SC  :  PWO.  Sea!  Beach  (  A 
NAVWPNSUPPCEN  C'ode  09  Crane  IN 
NCTC'  C'onst.  Elec.  School.  Port  Hucncmc,  CA 

NC'BC'  Code  10  Davisville.  Rl.  Code  15.  Port  Hueneme  C'A.  Code  155.  Port  Hueneme  C  A;  C Dde  15b.  Poit 
Huenemc.  C'a\;  C'ode  1571.  Port  Hueneme.  CA 
NC'BU  411  OIC.  Norfolk  VA 

NC  R  20.  C  ommander;  .XOth  Dct.  OIC.  Diego  Ciarcia  1 
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NMCB  74,  CO;  FIVE,  Operations  Dept 

NOAA  (Dr.  T.  N!c  Gumnessl  Rockville.  MI);  Library  Roek\ilie.  MI) 

NORDA  Code  410  Bay  St.  Louis,  MS;  Code  440  (Ocean  Rsch  Ofl)  Has  St  Louis  MS 
NRL  Code  58(K)  Washington.  DC;  C  ode  .>S4.^  (F  Rosenthal)  Washington.  DC.  C  ode  8441  (R.A.  Skop). 
Washington  DC' 

NROTC  J.W.  Stephenson.  I'C  .  Berkelev.  C  A 
NSD  SCE,  Suhic  Bay.  R  P. 

NTC  OICC,  CBl -401.  Great  Lakes  11 

NUCLEAR  RECiULATORV  COMMISSION  TC  Johnson.  Washington.  D( 

NUSC  Code  1.^1  New  London.  C^';  Code  332.  B-80  (J  Wilcox)  New  London.  CT;  Code  EA123  (R  S.  Munn). 

New  London  CT;  Code  TA131  (CL  De  la  Ou^).  New  London  CT 
OFFICE  SECRETARY  OF  DEFENSE  ASD  (MRA&L)  Code  (  SS  CX  Washington.  IX 
ONR  (Scientific  Dir)  Pasadena.  C'A;  Central  Regional  Offiee.  Bosiiin.  MA;  C\)de  481.  Ba\  Si  Louis.  MS; 

Code  485  (Silva)  Arlington,  VA;  Code  700F  Arlington  VA 
PHIBCB  1  P&E.  San  Diego.  CA;  1.  CO  San  Diego.  CA;  1.  C  SWC  D  Wellington.  San  Diego.  C  A 
PMTC  EOD  Mobile  L'nil.  Point  Mugu.  C'A.  Pat.  Counsel,  Point  .Mugu  CA 

PWC  CO  Norfolk.  VA;  CO.  (Code  Kb.  Oakland.  CA.  CX).  Great  Lakes  IL.  C’O.  Pearl  Harbor  HI;  Code  II). 
Great  Lakes.  IL;  Code  120.  Oakland  C'A;  C'ode  12()C.  (Library)  San  Diego.  C'A;  ('ode  128.  Guam;  Code 
154.  Great  Lakes,  IL;  Code  3l)C,  Norfolk.  V.A.  Code  .3()C.  San  Diego.  C'A;  Code  400.  Great  Lakes.  IL; 
Code  4(K).  Pearl  Harbor,  HI.  C'ode  4(K).  San  Diego.  C'A;  Code  420.  Cireat  Lakes.  IL;  Code  420.  Oakland. 
C'A;  Code  700.  Great  Lakes,  IL;  C  ode  XKI.  San  Diego,  C  A 
SUPANX  PWO.  Williamsburg  \'A 
UCT  ONE  QIC.  Norfolk.  VA;  OIC.  Port  Huenemc  CA 

L'S  DEPT  OF  INTERIOR  Bur  of  Land  Mngmnt  Code  583  (  L  F  Sullivan)  Washington.  DC 
C'S  GEOLOGICAL  SURVEY  Off.  Marine  Geoiogv.  Piteleki.  Resion  VA 
US  NAVAL  FORCES  Korea  (ENJ-P&O) 

USCG  (G-MP-3  USP/82)  Washington  Dc;  (Smith).  Washington.  DC' 

USCG  R&D  CENTER  CO  Groton,  C'T;  D.  Motherway.  Groton  CT;  Tech,  Dir.  Groton.  C'T 
USDA  Forest  Service  Reg  3  (R.  Brown)  Albuquerque.  NM 

USNA  Civil  Engr  Dept  (R.  Erchyl)  Annapolis  MD;  Ocean  Sys  Eng  Dept  (Dr.  Monney)  Annapolis,  MD;  PWD 
Engr  Div.  (C.  Bradford)  Anr)apolis  .MD 
WATER  &  POWER  RESOURCES  SERVICE  (Smoak)  Denver,  CO 


